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Abstract 

We report on microscopic numerical studies which support the chiral Lut- 
tinger liquid theory of the fractional Hall edge proposed by Wen. Our cal- 
culations are based in part on newly proposed and accurate many-body trial 
wavefunctions for the low-energy edge excitations of fractional incompressible 
states. 
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The quantum Hall effect (QHE) can occur only when a disorder-free two-dimensional 
electron system has an incompressibility, i.e., a discontinuity in the chemical potential, at 
a magnetic field dependent density, n*(B). The incompressibility implies a gap for charged 
and neutral excitations in the bulk of the system, but the magnetic field dependence of the 
density at which the gap occurs requires the existence of gapless excitations localized at the 
edge of the system. fl|] In equilibrium the current responsible for the orbital diamagnetism 
of the underlying two-dimensional electron system is carried at the edge and satisfies the 
thermodynamic identity, 

dl_ _ dn*(B) 

~dfl ~ ° dB ' ( ) 

In the edge-state picture, the quantum Hall effect follows from Eq. ( P when local equilibria 
are established on uncoupled edges. For the case of bulk incompressibilities at fractional 
Landau level filling factors (u), which are of many-body origin, Wen's chiral Luttinger liquid 
(CLL) theory |||J of the low-energy edge excitations predicts non-Fermi-liquid effects which 
have recently been of great interest. |3|-|5| The simplest version of this theory, and the one we 
test numerically here, applies to the edge of the incompressible ground states which occur 
at v = 1/m for m odd. [|5| 

The non-Fermi-liquid properties of one-dimensional fermion systems captured by the 
Luttinger model, |7j] arise from interactions between left-going and right-going particles in 
states close to the Fermi level. In a single-particle model, left-going and right-going states 
close to the Fermi energy in the QHE regime are localized on opposite edges of a sample with 
Hall bar geometry, and therefore interact weakly. However, Wen's theory predicts that in the 
fractional QHE regime, non-Fermi-liquid behavior arises without inter-edge interactions. We 
therefore consider a quantum Hall droplet (QHD) system for which electrons are confined 
to a finite area by a circularly symmetric external potential and for which a single circular 
edge exists. For such a system it follows directly from the microscopic fractional QHE 

physics !]] which gives rise to the chemical potential jump at v = 1/m, that the low-energy 
neutral edge excitations for any m are in one-to-one correspondence with those of a chiral 
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one-dimensional non-interacting fermion system which has single-particle states with only 
one sign of velocity. The non-interacting fermion system can in turn be mapped to a 
ID system of non-interacting chiral bosons; the bosons are the phonon modes of the edge. 
The same conclusion about the neutral excitations of a v — 1/m edge can be reached by 
using a hydrodynamic approach ]3|||] to derive a low-energy effective Hamiltonian expressed 
solely in terms of ID charge densities obtained by integrating the 2D charge density along 
the coordinate perpendicular to the edge. The Hamiltonian is quantized by invoking a 
commutation relation between Fourier components of the ID charge density: 

[p( q ),p(-q')] = ^S q , q/ . (2) 

For v = 1 Eq.( |2|) can be derived microscopically but for the fractional case it is assumed 
in order to satisfy Eq.( [1]) without altering the structure of the theory. This seemingly 
innocent introduction of the factor v on the right-hand-side of Eq.( |2|) is responsible for the 
non-Fermi-liquid behavior predictions [|3|-[§] of CLL theory. The predictions follow from the 
expression for the electron field operator, which plays a central role in the theory, in terms 
of boson operators. This expression, quoted below, is the simplest choice which results in a 
field operator that satisfies the correct commutation relation with the density operator. For 
the fractional case, it has not been possible to fully justify this expression on the basis of 
microscopic theory. This situation motivates the extensive numerical tests reported in this 
article. |H| Our findings are in complete agreement with the CLL theory. 



Our numerical calculations were performed for a QHD with a parabolic confinement 



potential. These model systems [|ITl|12f are of direct relevance to transport [13—15] and 



capacitance |16| measurements in quantum dots in strong magnetic fields. Here we are 



interested principally in the edge excitations of QHD's which are as large as possible in order 
to model the edge excitations of macroscopic incompressible states. In this model the single- 
particle state with angular momentum m has energy (m+ 1)Qq£ 2 where Qq is the frequency 
characterizing the parabolic confining potential and £ = (hc/eB) 1 ^ 2 is the magnetic length. 
For a QHD the angular momentum plays the role which would be played by the linear 



momentum in units of 2n/L in a Hall bar geometry. For this microscopic model system the 
incompressible ground states if>™(N) associated with the v — 1/m fractional QHE occur || 
for total angular momentum M = Mq{jti,N) = mN(N — l)/2 with N being the number 
of particles in the QHD. The m = 1 incompressible state is a single Slater determinant in 
which single-particle states with m = 0,1, ■ ■ ■ , N — 1 are occupied. 

In the language appropriate to the QHD the principal predictions of the CLL theory are 
the following: i) For M = Mo(m,N) + AM the spectrum has a low-energy branch with 
many-particle eigenenergies given by E = Eq + Yli^i^i where J2i^ n i — AM and Eq is the 
energy of the Mo(m,N) incompressible state. This property is expected to be accurately 
satisfied for AM < N 1 ^ 2 , since the excitations are then well localized at the edge. Here rii 
are the non-negative integers which give the occupation numbers for the bosonic edge-wave 
angular momentum / and energy e\. ii) The electron creation operator is given by 



a] is a boson creation operator, z is a constant which is not fixed in the CLL theory, 
4>-(0) = [</>+(#)] t, and for each particle number the boson operators act on the bosonic 
quantum numbers of the edge waves. 

We first discuss our tests of the CLL theory predictions for the bosonic nature of the 
excitation spectrum of the edge of the QHD. Substantial arguments can be advanced in favor 
of this aspect of the CLL theory predictions from microscopic theory. In the case of the 
hard-core model of electron-electron interactions, the low-energy portion of the spectrum has 
no contributions from electron-electron interactions, the bosonization follows from analytic 
arguments and e/ = IQli 2 . More generally, qualitative aspects of the bosonization of 
the low-energy portion of the spectrum follow from the adiabatic evolution of the spectrum 
with changing model interactions. For the QHD the / = 1 single-boson state corresponds 



^(0) = ^exp[0 + (0)]e X p[-0_(0)] 



(3) 



where 




(4) 
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microscopically to an excitation of the center-of-mass of all electrons from M cm = to 
M cm = 1. Since the interaction energy of the electrons is independent of the center-of-mass 
state [|TT| , |T8f it follows e\ = Qq£ 2 , independent of electron-electron interactions. For I 7^ 1 ei 



is interaction-dependent; to test the bosonization for the physically realistic Coulomb inter- 
action model we have determined the spectrum of finite-size QHD's by exact diagonalization 
of the many-particle Hamiltonian neglecting confinement. For parabolically confined QHD's 
the eigenstates are unchanged and the subspace spectrum at total angular momentum M 
shifts rigidly by MQq£ 2 when the confinement is introduced. 

Figure [I] shows the spectra of the QHD for Qq = as a function of M close to 
M = M (3,N) and M = M (1,N). Note that the interaction energy decreases as M 
increases because of the decrease in average two-dimensional electron density; in the QHD 
case interactions lower the boson energies because the charge spreads out in the direction 
perpendicular to the edge when the angular momentum is changed. Fig.( p]) shows that 
the bosonization law for the neutral edge-wave excitation spectrum, which is exact for the 
hard-core model, is still closely obeyed for the physically realistic Coulomb interactions. The 
bosonization is even more robust than would be expected a priori since it appears to hold 
even where e; oc / fails. 

More extensive analyses of the spectrum are possible and have been completed previ- 
ously only for the edge excitations of a v = 1 QHD where exact diagonalization techniques 



can be applied for much larger numbers of electrons. [[L9|] A set of microscopic operators 



have been proposed by Stone |19| to generate the edge-wave spectrum in the v = 1 case: 



&AM = E m (( m + ^^)V m ') 1 ^ 2c m+AAf c m where and c m are the electron creation and 



annihilation operators. Recently, Oaknin et al. [|20[ identified a modified set of operators, 
4m = Em( m '/( m + / ^^)0 1 ^ c L+AM c m which, when acting upon the v = 1 QHD, do not 
alter the center-of-mass state of the electrons. It was demonstrated by Oaknin et al. that 
their operators generate the single-boson excitations of the v = 1 QHD more accurately 
than those of Stone. [^] More generally we find by comparing with exact eigenstates, that 
for a given AM the excitations with large bosonic occupation numbers are more accurately 



generated by the S* operators than by the operators, while the comparison is reversed 
for small bosonic occupation numbers. In the limit AM < N l l 2 the states generated by the 
two sets of operators become equivalent. 

In order to access large N for a fractional QHD, following Ref. || we propose a set of 
microscopic many-body wavefunctions for low-lying excited states of fractional QHD's: 

D m - x n {j}r (5) 
i 

where X)/ l n i — AM and D is the Vandermonde determinant f| Jastrow factor which relates 
different Laughlin QHD states [ty™(N) = D 2 ^™~ 2 (N)}. (As discussed above for some states 
greater accuracy at finite N can be achieved on substituting Jt by S^). These trial wave 
functions have the following properties: (i) they have the appropriate value of the angular 
momentum; (ii) for AM = they reduce to Laughlin's approximation for the incompressible 
state many-body wavefunction; (iii) they are exact eigenstates with zero interaction energy 
in the case of the hard-core model Hamiltonian; and (iv) they are eigenstates of center-of- 
mass angular momentum with eigenvalue for AM ^ 1 and eigenvalue 1 for AM = 1. 
If the set S' is used instead Jt the two final properties still hold but only in the limit 
AM < N 1 / 2 . Table I shows the absolute value of the overlaps between the trial single-boson 
state wavefunctions and the exact ones for the v — 1/3 QHD for up to six particles. We see 
that the excited state wavefunctions (AM > 2) tend to be, if anything, more accurate than 
the Laughlin trial wavefunction for the ground state (AM = 0). (The bracketed overlap 
values were calculated using the Stone operators). It is important to realize that the direct 
application of the Oaknin et al. operators to the v = 1/3 Laughlin state would not generate 
identical states since [D" 1 ^ 1 , j\ M \ ^ 0; the resulting states would lack properties iii) and iv) 
from the list above and they are very poor approximations to the true states. On the other 
hand [D m -\Si M ] = 0. 

We turn finally to our numerical tests of the less fully justified CLL theory expression for 
the the electron field operator. We report here only results for the squared matrix elements 
connecting the ground state of the iV electron system with the ground state and low-lying 
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bosonic states of the iV + 1 particle system at v = 1/3, \(^ 3 {ni} (N + l)\^(9)\^l(N))\ 2 . 
The values predicted for these squared matrix elements by the CLL theory can easily be 
computed from Eq. (|3]). The predicted matrix elements for AM = 0, 1,2,3,4 are listed in 
Table II and compared with those calculated microscopically |H[ for N = 6,7, and 8. For 
each case the matrix elements have been normalized to the ground-state-to-ground-state 
matrix element, z. |22| Since the angular momentum difference between the N and N + 1 
particle states is M (3, iV + 1) + AM - M (3, N) = 3N + AM only the part of the electron 
creation operator proportional to c\ N+AM contributes to the microscopic matrix element. 
The agreement between the CLL results for these matrix elements and the microscopic 
calculations is excellent. It appears from our extrapolation that the CLL theory becomes 
exact for N — > oo. The non- Fermi-liquid power law properties predicted by the CLL theory, 
depend both on the the predictions for these matrix elements and on having dispersionless 
propagation of edge modes, i.e., ei = cl. In that case the total spectral weight at an energy 
cl above the chemical potential is proportional to the sum of the squared matrix elements 
for AM = I. It follows from Eq. (^) that this sum equals, 

(AM + m-1)! 
Aam = Z (AM)!(m-l)r (6) 

For AM ^> (m — 1) the right-hand side of Eq.( [|) approaches (AM)™" 1 / (m — 1)!. This gives 
the power law dependence of the spectral weight at low energies which enters into various 
physical properties. 

In conclusion, our microscopic numerical studies strongly support the chiral Luttinger 
liquid theory of the fractional Hall edge proposed by Wen and add to the motivation for 
experiments which can probe the low-energy excitations of incompressible fractional Hall 
states. We thank Luis Brey, Matthew Fisher, Steve Girvin, Rudolf Haussman, Charles Kane, 
Kyungsun Moon, Michael Geller, Jacob H. Oaknin, Carlos Tejedor, and Ulrich Zuelicke for 
informative and stimulating discussions. This work has been supported by the National 
Science Foundation under grant DMR-9416902. J.J.P acknowledges support from NATO 
postdoctoral research fellowship. 
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TABLES 

TABLE I. Absolute value of the overlaps between the exact single-boson state (see text) and 
the trial states D 2 J am \^q(N)) for up to six particles (in brackets using S AM ). The high overlaps 
confirm the validity of the mapping (through D 2 ) between the single-boson state excitations of a 
v = 1 QHD and those of a v = 1/3 QHD. 
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TABLE II. Microscopic and CLL theory spectral weights for the v = 1/3 QHD state in units 
of the ground-state-to-ground-state matrix element. Results are shown for N = 6, 7, 8 and are 
extrapolated to N — > oo. The bracketed CLL theory results are for v = 1 where the squared 
matrix elements sum to one for each AM. For v = 1/3 each matrix element is increased by a 
factor of 3 fc where k = J2i n i- The increase in the spectral weights moving away the Fermi level is 
due to the increasing boson occupation numbers. 
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FIGURES 

FIG. 1. Low-lying excitation spectra of a 6-particle QHD neglecting confinement energies, (a) 
v = 1 The incompressible ground state occurs at M=15. The solid dots show the interaction 
energies of the single-boson edge states, e\. [e\ = 0) The highlighted asterisk at M=19 shows the 
energy of the state with ri2 = 2. Its energy differs from 2e2 by 3.3%. The sequence of states 
which appear along horizontal lines in this figures correspond microscopically to increasing values 
of the center-of-mass angular momentum and, in the boson picture, to increasing values of rt\. (b) 
v = 1/3 The incompressible ground state occurs at M=45. The similarity with the u = 1 spectrum 
in (a) is clear despite the appearance for AM > 3 of other states representing bulk excitations. 
These states are expected to move to relatively higher energies for large N. If higher LL's were 
included in the calculation, the same incursion of bulk excitations for values of AM near N could 
occur in the v = 1 case at small enough values of hco c . In the v = 1/3 case the energy of the 
highlighted ri2 = 2 state differs from 2e2 by 0.3%. 
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